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Abstract. We derive the expression for the snap parameter in f(R) gravity. We use 
the Palatini variational principle to obtain the field equations and regard the Einstein 
conformal frame as physical. We predict the present-day value of the snap parameter 
for the particular case f(R) — R— const/ R, which is the simplest f(R) model explaining 
the current acceleration of the universe. 
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1. Introduction 

f(R) gravity models, in which the gravitational Lagrangian is a function of the curvature 
scalar R [fl |2], have recently attracted a lot of interest. They explain how the current 
cosmic acceleration originates from the addition of a term R~ l to the Einstein-Hilbert 
Lagrangian R [3j. As in general relativity, f{R) gravity theories obtain the field 
equations by varying the total action for both the field and matter, and equaling 
this variation to zero. Here we use the metric-affine (Palatini) variational principle, 
according to which the metric and the affine connection V ^ v are considered as 
geometrically independent quantities, and the action is varied with respect to both of 
them [H El [7]. The standard approach is the metric (Einstein-Hilbert) variational 
principle, according to which the action is varied with respect to the metric tensor, and 
the affine connection coefficients are the Christoffel symbols of g^ v [8]. Both approaches 
give the same result only if we use the standard Einstein-Hilbert action [3]. The field 
equations in the Palatini formalism are second-order differential equations, while for 
metric theories they are fourth-order [9]. Another remarkable property of the metric- 
affine approach is that the field equations in vacuum reduce to the standard Einstein 
equations of general relativity with a cosmological constant [91 [TO] . 

One can show that f(R) theories of gravitation are conformally equivalent to the 
Einstein theory of the gravitational field interacting with additional matter fields, if 
the action for matter does not depend on the connection [7J O [11] . This can be done 
by means of a Legendre transformation, which replaces an f(R) Lagrangian with a 
Helmholtz Lagrangian [3 [12]. For f(R) gravity, the scalar degree of freedom due to the 
occurrence of nonlinear second-order terms in the Lagrangian is transformed into an 
auxiliary scalar field <fi [11]. The set of variables (g^u, 0) is commonly called the Jordan 
conformal frame. In the Jordan frame, the connection is metric incompatible unless 
f(R) = R. The compatibility can be restored by a certain conformal transformation 
of the metric: g^ v — > h^ v = f'(R)g IMU [13J. The new set (/i M „, (p) is called the Einstein 
conformal frame, and we will regard the metric in this frame as physical. Although 
both frames are equivalent mathematically, they are not equivalent physically [H] , and 
the interpretation of cosmological observations can drastically change depending on the 
adopted frame [T5] . 

f(R) gravity models in both the metric and metric-affine formalism have been 
compared with cosmological observations by several authors [HH [TFl [TEl [191 ED] • The 
problem of viability of these models at smaller scales, namely their compatibility 
with solar system observations, is a subject of recent debate ([21] and references 
therein). There are also limits on these models arising from particle physics [22] , 
matter instability [23] and violation of the equivalence principle [14, 24J. Current SNIa 
observations provide the data on the time evolution of the deceleration parameter q in 
the form of q = q{z), where z is the redshift [25 j . The extraction of the information 
from these data depends, however, on assumed parametrization of q(z) [26]. For small 
values of z such a dependence can be linear, q(z) = q + q 1 z [25J, but its validity should 
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fail at z ~ 1. A convenient method to describe models close to ACDM is based on the 
cosmic jerk parameter j, a dimensionless third derivative of the scale factor with respect 
to the cosmic time [27J EE] • A deceleration-to-acceleration transition occurs for models 
with a positive value of jo an d negative qo.S The flat ACDM models have a constant 
jerk .7 = 1. 

Recently, we showed [29J that the predictions for the current value of the jerk 
parameter for the particular case f(R) = R—f^, which is the simplest way of introducing 
the cosmological term in f(R) gravity [21 ED], agree with the SNLS SNIa |30J and x- 
ray galaxy cluster distance [26] data, but do not with the SNIa gold sample data [25J. 
Moreover, the predicted value of the deceleration parameter in this model agrees with all 
three data sets [29J. Therefore f(R) models based on the Palatini variational principle 

2 

and formulated in the Einstein frame, including the case f(R) = R — f^, provide 
a possible explanation of the current cosmic acceleration. More constraints on these 
models are likely to come from nondimensional parameters containing higher derivatives 
of the scale factor, such as the snap parameter s = [28] . 

In this work, we find the general expression for the snap parameter in f(R) gravity, 
assuming that the universe is flat, homogeneous, isotropic and pressureless. We use the 
field equations derived from the Palatini variational principle. We assume that matter 
is minimally coupled to the metric tensor in the Jordan frame, and then transform to 
the Einstein frame which we consider physical [6] and in which this coupling is non- 
minimal Since the question of which frame is physical and in which frame the 
matter-metric coupling is minimal should be ultimately answered by experiment or 
observation, the presented model should be treated as phenomenological. Anticipating 
cosmological measurements, we predict the current value of the snap parameter for the 
case f(R) = R - fj. 

2. The field equations in f(R) gravity 

The action for f(R) gravity in the original (Jordan) frame with the metric g^ u is given 
by [6] 

Sj = ~ I #x[f^f(R)} + S m (9^, # (1) 

Here, y/—gf(R) is a Lagrangian density that depends on the curvature scalar: R = 
R^viT p a )g iiv , S m is the action for matter represented symbolically by if) and independent 
of the connection, and k = Tildes indicate quantities calculated in the Jordan 

frame. 

Varying the action Sj with respect to g^ u yields the field equations: 

f\R)R, p - \f{R)g„ = kT, u , (2) 

where the dynamical energy-momentum tensor of matter, T^, is generated by the 
Jordan metric tensor [6] : 

6Sm = YcI d4x ^~ Zg T ^ 5 9^> ( 3 ) 
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and the prime denotes the derivative of a function with respect to its variable. From 
vanishing of the variation of Sj with respect to the connection V? it follows that 
this connection is given by the Christoffel symbols of the conformally transformed 
metric pH [13] 

9,u = f\R)g^ (4) 

The metric defines the Einstein frame, in which the connection is metric compatible. 
The action (CQ) is dynamically equivalent to the following Helmholtz action [6ll9j[TT]: 

1 r + 



= P- (6) 

R=<f,(p) 



Sh = -— d i x^-~g[f(<f ) (p))+p(R-<f ) (p))] + S m (~9,u^), (5) 
where p is a new scalar variable. The function (p(p) is determined by 

df(R) 
dR 

From equations (HI) and (EJ) it follows that 

<P = Rf'(<P), (7) 

where R = R lll/ (Tp a )g IJiU is the Riemannian curvature scalar of the metric g^ v . 

In the Einstein frame, the action (jSJ) becomes the standard Einstein-Hilbert action 
of general relativity with an additional scalar field: 

S E = ~ I afx^-giR-p^m+P^f^m + Sm^g^)- (8) 



2kc . 

Choosing (which is equal to the curvature scalar R in the Jordan frame) as the scalar 
variable leads to 

Se = -^J d i x^[R-2V{<j>)] + S m {[f\4>)]' 1 g^l (9) 
where V(<f>) is the effective potential: 

w Arm 2 ' 

Varying the action ([9]) with respect to g^ yields the equations of the gravitational 
field in the Einstein frame [5j [6] : 
1 kT 

R^v - ^Rg^v = - Viflguv, (li) 

while the variation with respect to <fi reproduces (|7J). Equations (171) and (ITTil give 

<j>f'(<j>)-2f(<f>) = KTf'(<j>), (12) 

from which we obtain = 4>{T). Substituting <fi into the field equations (fTTj) leads to a 
relation between the Ricci tensor and the energy-momentum tensor, 

R,xu - ^Rg^ = K r (T)T^ + A(T)g^, (13) 

with a running gravitational coupling K r (T) = k[/'((/)(T))] _1 and a variable cosmological 
term A(T) = -V(<f>(T)): 

m = m=*fM. (14) 
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= + HtiFTr. (17) 



Such a relation is in general nonlinear and depends on the form of the function f(R). 
The Bianchi identity applied to equation ffTTl) gives 

V 1 (0) V f ( )+ 2 k[/'(0)]' > (15) 

This relation means that the energy-momentum tensor in the Einstein frame is not 
covariantly conserved, unless f(R) = R or T = |20j. We can write the field 
equation (TTTT) as 

R^-\r9^ = <T^ + T%), (16) 

where T™ v = T^. This defines the dark energy-momentum tensor, 

AM 1- W ), 
« ^ f(0) 

From equation ffTBT) it follows that matter and dark energy form together a system that 
has a conserved 4-momentum. Consequently, in the Palatini f(R) gravity formulated in 
the Einstein frame, matter and dark energy interact ([31] and references therein). This 
interaction may be responsible for the observed large entropy of the universe. 

We assumed that matter is minimally coupled to the metric tensor in the Jordan 
frame. Then we transformed to the Einstein frame, in which this coupling becomes 
non-minimal, and assumed that this frame is physical, motivated by the fact that 
the connection is metric compatible in this frame. Such a construction is completely 
phenomenological. However, if we consider the Einstein frame from the very beginning 
and define the energy-momentum tensor generated by the metric tensor g^ v as the 
true energy-momentum tensor for matter (minimal coupling in the Einstein frame), the 
resulting action, instead of equation (jHJ), is 



S E = / d 4 x^g~[R - 2V {</>)] + S m {g^ # 



Consequently, the field <fi does not couple to anything and we arrive at general relativity 
with the cosmo logical constant, which is not interesting from a modified gravity 
perspective [HE]. 

3. The snap parameter in f(R) gravity 

The snap parameter in cosmology is defined as [28] 

- < i9 > 

where a is the cosmic scale factor, H is the Hubble parameter, and the dot denotes 
differentiation with respect to the cosmic time. This parameter appears in the fourth- 
order term of the Taylor expansion of the scale factor around ao: 

^ = 1 + H (t - t ) - \q*Hl{t - t ) 2 + \j,Hl{t - t ) 3 

GO AO 

+ ^soH 4 (t - t ) 4 + 0[(t - t ) 5 ], (20) 
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where the subscript denotes the present-day value. We can rewrite equation (|TU|) as 

s=jj-j(2 + 3q), (21) 

where q is the deceleration parameter and j is the jerk parameter. For the flat KCDM 
model s = — (2 + 3g) since j = 1 [261 [29], an d the departure of the quantity ds/dq from 
—3 measures how the evolution of the universe deviates from the ACDM dynamics. 

From the gravitational field equations ( TTTT) applied to a flat Robertson-Walker 
universe filled with dust we can derive the 0-dependence of the Hubble parameter [6] 



the deceleration parameter [20 



_ 20f(0)-3/(0) 

qW ~ 4>m - 3/(0) ' (23J 

and the jerk parameter [2H] 

m = wr + io0 3 r 3 r - 750 2 /' 2 //" - 3 + is/ 2 /' 2 - i62/ 3 r 

+ 1890/ 2 /'/"] x W - 3/) 2 (2f 2 + 0/7" - 6//")]- 1 . (24) 

Th prime denotes the differentiation with respect to 0. We also have the expression for 
the time dependence of [6]: 

• Vec^f - vww - 3/ 

9 2f>* + <f>f>f»-6ff» • 1 } 

For the 0-derivative of the jerk parameter we obtain a quite complicated expression: 

f = Uf - 3/)(2f 2 + 0/7" - 6/f )(3O0 3 f 2 f 2 + 1O0 3 /' 3 f' 

- 15O0 2 //7" 2 - 370 2 f 3 f - 750 2 //' 2 f' - 80f 4 + 24/f 3 + 1890/ 2 f 2 
+ 1890/ 2 /'/"' + 1920/f 2 /" - 162/ 3 f - 267 fff") 

- (20 2 f 4 + 1O0 3 /' 3 f - 750 2 f 2 /f - 120/f 3 + 18/ 2 f 2 + 1890/ 2 ff 

- 162/ 3 /") x (30 2 /7" 2 - 150//" 2 - 8f 3 + 27//7" - 0f 2 f + 2 f 2 f 

- 90///"' + 18/ 2 /'")] x [(0f - 3/) 3 (2f 2 + 0/7" - 6//") 2 ]- 1 - (26) 
Combining equations fl2~lT[25l) and using j = 0j'(0) lead to 

= ., 6fW - 2f) _ . 80/' -15/ 
S 3 (2f 2 + 0f/" - 6ff") 3 0/' - 3/ ' 1 ] 

Putting here j from equation(T54|) and j' from equation ff26l) gives the final expression for 
the snap parameter in f(R) gravity as a function of 0, /(0), /'(0), /"(0), and /'"(0), 
which we do not write explicitly. 

We now examine the case /(i?) = R — f^, where a is a constant, which is a possible 
explanation of the current cosmic acceleration [3]. In this model, the present-day value 
of is 0o = (-1.05 ± 0.01)a, where a = (7.35tJ;l?) x lO" 52 ™" 2 [20]. We do not 
need to know the exact value of a since it does not affect nondimensional cosmological 
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parameters. Substituting O into equations (1211) . (125]) and (12"T|) gives the present-day 
value of the cosmic snap parameter^] 

so = -0.22±8;fJ. (28) 

2 

In the /(-R) = R — fx model, the deceleration-to-acceleration transition occurred at 
4>t = —\Jh/?>a [20]. Consequently, we find the snap parameter at this moment: 

s t = -2.68. (29) 

This value shows that the snap parameter in f(R) gravity changes significantly between 
the deceleration-to-acceleration transition and now, which is clear from equation (|21j) 
and the fact that the deceleration parameter changes in this period of time from to the 
predicted value q = — 0. 67io!o3 [2D]- For the flat AC DM model, the snap parameter 
increases from s = —7/2 for the matter epoch, through s = —2 at this transition, to 
the asymptotic de Sitter value s — 1, indicating the difference between the f(R) and 
KCDM predictions for s t . 

Lastly, we show the role of the energy conditions [32] in metric-affine f(R) gravity 
models. For a pressureless universe, these conditions reduce to the inequality 

e > 0, (30) 

where e = T is the energy density of matter. From equations (Tl2l) . ( J22l) and ( J23l) we 
obtain 

e = **£±MM. (31 , 

Since q > — 1 (it reaches —1 asymptotically [2D]) and /'(</>) > (this condition assures 
that the conformal transformation from the Jordan to the Einstein frame does not 
change the signature of the metric tensor), formula (1301) is satisfied. Therefore, the 
energy conditions do not impose additional constraints on Palatini f(R) gravity models. 
For metric f(R) models, these conditions lead to constraints containing the jerk and 
snap parameter [33]F 



4. Summary 

We derived the expression for the cosmic snap parameter in f(R) gravity formulated 
in the Einstein conformal frame. We used the Palatini variational principle to obtain 
the field equations and apply them to a flat, homogeneous, and isotropic universe filled 
with dust. We considered the particular case f[R) = R — f^, which is the simplest f(R) 
model explaining the current cosmic acceleration, and for which the predicted present- 
day values of the deceleration and jerk parameters are quite consistent with cosmological 
data. For the present-day value of the snap parameter, we predict sq = — 0.22i ;fg. 

| The predicted value for the current cosmic jerk parameter found in [29] is jo = I.OI^q'^i- Here, we 
recalculated this value and obtained jo — 1.01 ± 0.01, which differs from the former by the precision 
errors. This correction does not change the conclusions of [29] . 

§ The jerk and snap parameters do not appear in equation (|3ip and the energy conditions because the 
field equations are second order, as opposed to the fourth-order field equations in the metric formalism. 
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Expanding the scale factor to fourth order with respect to time (equation (I2UD ) is 
physically meaningful, since cosmological data already allow to measure the third-order 
term: the jerk parameter [26]. The snap parameter may be important for observations 
involving redshifts z ~ 1 and higher, where expansion of cosmological quantities in 
powers of z cannot be limited only to linear and quadratic terms. Therefore this paper 
is not only a formal mathematical exercise, but also provides physically measureable 
constraints on Palatini f(R) gravity. 

A cosmological sequence of matter dominance, deceleration-to-acceleration 
transition and acceleration era may always emerge as cosmological solutions of f{R) 
gravity [17J. We showed in [31] that, in Palatini f(R) gravity, the deviation of the 
growth of the cosmic scale factor in the matter era from the standard law a{t) ~ t 2 ^ 
is small, which is consistent with WMAP cosmological data |34j . On the other hand, 
metric f{R) gravity models with a power of R dominant at large or small R yield the law 
a(t) ~ t 1 / 2 which is ruled out by cosmological observations [35]. Therefore f(R) gravity 
in the Palatini variational formalism is a viable theory of gravitation that explains the 
current cosmic acceleration. 

It is possible to find an f(R) action without cosmological constant which exactly 
reproduces the behavior of the Einstein-Hilbert action with cosmological constant, 
i.e. the expansion history of the universe does not uniquely determine the form 
of the gravitational action (T8J [19]. Moreover, the background expansion alone 
cannot distinguish between different choices of f(R) and one must study cosmological 
perturbations in order to determine which choice is physical [191 EE]. Measurements 
of higher derivatives of the scale factor, including the snap parameter, will probably 
constitute a robust part of these studies. 
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